Chirped chiral solitons in nonlinear Schrodinger equation with self-steepening and 

self-frequency shift 
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We find exact solutions to nonlinear Schrodinger equation in the presence of self-steepening and 
self- frequency shift. These include periodic solutions and localized solutions of dark-bright type 
which can be chiral, and chirality being controlled by sign of self steepening term. A new form 
of self phase modulation, which can be tuned by higher order nonlinearities as also by the initial 
conditions, distinct from nonlinear Schrodinger equation, characterizes these solutions. In certain 
nontrivial parameter domain solutions are found to satisfy linear Schrodinger equation, indicating 
possiblity of linear superposition in this nonlinear system. Dark and bright solitons exist in both 
anomalous and normal dispersion regimes and a duality between dark-bright type of solution and 
kinematic-higher order chirping is also seen. Localized kink solutions similar to NLSE solitons, but 
with very different self phase modulation, are identified. 
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Nonlinear Schrodinger equation (NLSE): 
itpx + aiiptt + a2|V'lV = 0, 



(1) 



governs the dynamics of picosecond pulse propagation in 
optical fibers [I|], where ai is group velocity dispersion 
(GVD) parameter and 02 specifies the strength of Kerr 
nonhnearity. As predicted by Hasegawa and Tappert J2| 
and experimentally observed by MoUenauer et al. [3|, 
this system supports stable soliton solutions owing their 
existence to complete integrability [3]. With the advent 
of high intensity laser beams, it has become possible to 
generate optical pulses with width of the order of 10 fem- 
toseconds. The higher order effects like third order dis- 
persion, self steepening of pulse due to dependence of the 
slowly varying part of nonlinear polarization on time and 
self frequency shift arising from delayed Raman response 
become important while studying propagation of these 
pulses. Inorder to account for them, Kodama 5] and 
Kodama & Hasegawa H proposed higher order nonlin- 
ear Schrodinger equation (HNLSE) as a generalization of 
NLSE: 



(2) 



where a third order dispersion with coefficient 03, self 
steepening term with coefficient 04 and self frequency 
shift effect with coefficient 05 have been added. This 
model, unlike NLSE, is not integrable in general. A 
few integrable cases have been identified: (i) Sasa- 
Satsuma case (a3:a4:(a4-|-a5) = 1:6:3) [ij, (ii) Hirotacase 
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(a3:a4:(a4 + as) = 1:6:0) ^ and (iii) derivative NLSE of 
type I and type II 0]. Many restrictive special solutions 
of bright and dark type have been obtained [l^, [HI, • 
The effect of third order dispersion is significant for 
femtosecond pulses when GVD is close to zero. It is neg- 
ligible for optical pulses whose width is of the order of 
100 femtoseconds or more, having power of the order of 
1 Watt and GVD far away from zero. However, in this 
case self steepening, as well as self frequency shift terms 
are still dominant and should be retained. The effects of 
these higher order terms on pulse propagation have been 
extensively studied numerically |l|, Il3l | , and some special 
solutions to this system are also known [T^ . 

In this letter, we report existence of a new class of lo- 
calized, as well as periodic solutions for NLSE in the pres- 
ence of self steepening and self frequency shift. The lo- 
calized solutions include dark, bright and kink type soli- 
tons. These complex solitons are generally chirped and 
show a new form of self phase modulation (SPM). Unlike 
NLSE, where chirping is of kinematic origin (controlled 
by initial conditions), the chirping in present case has 
both kinematic as well dynamic origin. The former varies 
as reciprocal of intensity, whereas the latter is directly 
proportional to intensity and depends upon higher order 
nonlinearities only. It is evident that kinematic chirping 
plays a significant role for dark soliton, where as the dy- 
namical chirping would be important for bright soliton. 
We found that for certain special values of parameters the 
intensity profile of these solitons are similar to NLSE soli- 
tons, while retaining the nontrivial phase structure due 
to higher order nonlinear terms. It is also seen that, in 
the same regime, kink solitons exist which otherwise are 
not allowed in this system. For certain parameter values, 
this model mimics NLSE, however the presence of higher 
order nonlinearities crucially affects the dynamics. And 
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as a consequences one sees that fundamental bright soh- 
ton of NLSE is no longer a valid solution. Further, trivial 
phase dark solitons exist in anomalous dispersion regime, 
are found to be chiral, with the direction of propagation 
set by self steepening term. Chirped dark and bright 
solitons are found to exist in normal as well as anoma- 
lous dispersion regimes. Chirped dark solitons in anoma- 
lous regime, whereas chirped bright solitons in normal 
dispersion regime exhibit chirality, which is controlled 
by self steepening term. Observation of chiral solitons 
is one of the main results of this letter. Very interest- 
ingly, it is seen that for some nontrivial choice of param- 
eters, this system behaves like free particle Schrodinger 
equation with appropriate constants, making this system 
amenable to linear superposition, which is otherwise for- 
bidden in this nonlinear system. This system is found to 
be Painleve integrable thereby establishing existence of 
regular solutions. 

Modulo a trivial kinematic phase, the complex enve- 
lope travelling wave solutions can be generally repre- 
sented as: 



ipix,t) = p(Oe 



(3) 



where ^ = a{t — ux) is the travelling coordinate, and p 
and X are real functions of ^. Here, a is scale parameter 
and u = 1/v with v being the group velocity of the wave 
packet. The ansatz solution leads to the compatibility 
conditions: 



— aup' + 2a^aix' p' + a^aix" p 
^iaa^p^p' + 2aa5p^p' = 0, (4) 
aux' p -\- c?a\p" — a^aix'^ p -\- a2P^ — aa^x^ p^ = 0. (5) 

Equation (U) can be exactly integrated to yield: 



X 



2aai 



aaip^ 



(3a4 + 205) 2 



(6) 



where c is to be determined by initial conditions. It is to 
be noted that the phase has a nontrivial form and has two 
intensity dependent chirping terms, apart from kinematic 
first term. As is evident the second term is of kinematic 
origin and is common to Schrodinger equation as well. 
The last term is due to higher nonlinearities and leads 
to chirping that is exactly inverse to that of the former. 
This is a novel form of self phase modulation which is 
controlled by interaction. The amplitude equation ([5]) 
reduces to: 



hp" + O2P + Osp^ + I 



p^ 



(7) 



with Ol = a'ai', 62 = (" -ca.+2ca,) ^ ^ {2a^a^-ua^) 

and 6I4 = (4°4-i)(3a4-H2a5) ^ ^^^^ ^j^^^ ^j^^ uoutrivial 
contribution from higher order nonlinear terms is through 
04, which is zero for a4 — 1/4 or 04 : 05 = —2 : 3 
(assuming 04 7^ and 05 7^ 0), and as we shall soon 



show, this results in interesting physical consequences. 
In the case when 04 : 05 = — 2 : 3, both the intensity 
as well as phase will not be having any new features due 
to higher order nonlinearities and the solutions will ex- 
actly resemble NLSE solutions. However unlike NLSE 
both dark and bright solitons exists, in both normal and 
anomalous dispersion regime. When c = 0, existence of 
fundamental bright solitons with p = ^sech(^) is forbid- 
den, and only dark solitons with p =^ ^tanh(^) exist, 
with a = T^^. Furthermore, these dark solitons in 

zai ^ ' 

anomalous dispersion regime, i.e., ai > 0, respect the 
inequality a^u > 2|ai||a2|, which restrict them to travel 
only along one direction, given by sign of a^. This is 
an example of chiral soliton which is absent in NLSE. 
When 0, both bright soliton and dark soliton exists, 
and satisfy 2aia2 > 1*04 and 2aia2 < ua4 respectively, 
showing that both have mutally exclusive velocity space. 
As a consequence, in anomalous dispersion regime, dark 
solitons obey the inequality 2|ai||a2| < ua^ hence are 
chiral, whereas in normal dispersion regime bright soli- 
tons satisfy — 2|ai||a2| > ua^ and hence are also chiral. 
Notice that the directionality of these solitons is due to 
the presence of higher order terms, the sign of decides 
the direction in which solitons are allowed to propagate. 
For a4 = 1/4, the intensity profile will be the same as 
NLSE, whereas phase will still show nontrivial chirping. 
In this case also both dark and bright solitons can be 
chiral, and can exist in normal and anomalous disperison 
regimes, which is in sharp contrast to NLSE. 

It is very intriguing to see that when 04 = — and 
u — , equation ^ combined with equation ([6]) re- 

duces to free particle Schrodinger equation in ■(/;. So for 
this choice of parameters, in the presence of both Kerr 
and higher order nonlinearities, the effective evolution 
equation for is linear, and one would expect to see phe- 
nomena like interference, which is forbidden otherwise in 
this system. 

Equation ([7]) can be cast into a convenient form using 
p^^/a: 
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(8) 



where k is constant fixed by initial conditions. Solutions 
for this equation, with 64 ^ 0, can be found by conformal 
Mobius transformation: 



A + Bf 
C + Df 



(9) 



which for some suitable A, B, C and D, connects cr to 
elliptic function /. These elliptic functions, as is known, 
are generalization of trigonometric and hyperbolic func- 
tions and appear in solutions of many nonlinear equa- 
tions. 

Considering the importance of localized solutions, we 
set /(^) = sech{£_), and look for allowed values of A, B, 
C and D for which ^ is solution of equation ([H]). The 
consistency conditions leads to: A = 8^4 A — 36*3, B = 
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FIG. 1: Intensity profile of few solutions: i) Dark Soliton (in 
red color) for a'l = 1.6001, a2 = -2.6885, as = -2.8302, 
ai = 0.30814, as = 0.76604, u = 4.1185, c = -3.1186, k = 
-77.965, A = 35.36, a = 4.7421, B = -11.912, C = 4.5702 
and D — 17.855; ii) Kink Soliton (in blue color) ai = —11.197, 
a2 = 44.778, 03 = 6.219, 04 = 19.066, as = 37.301, u = 
884.36, c = -13810, fc = 4360.2, A = 0.94014, a = 8.6548, 
B = 0.06745 and D = -0.97921. 



864,8 -W3D, C = 804 and D 
a are given by: (10246l/)i3- 



86*4 D where A, B, D and 
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189, 



= ±^ 

V54e3" -27933-969294-6494^^2 

-9932-329294 -64A942 



and 
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Since the exact closed form solution is known, a simple 
maxima-minima analysis is sufficient to distinguish pa- 
rameter regimes supporting dark/bright solitary waves 
p^ . In this case, when AD > BC one gets a bright 
soliton, whereas if AD < BC then dark soliton exists. 
Figure (P) depicts intensity profile of a typical dark soli- 
ton. It is interesting to note that for dark solitons, 
from equation ([6]), the kinematic chirping is dominant 
at the centre of the pulse whereas the higher order chirp- 
ing is dominant away from the center (see figure 
However, exactly the opposite is true for bright solitons 
where the center is dominated by higher order chirp- 
ing and kinematic chirping is important away from the 
center. This shows that there is duality between dark 
soliton-kinematic chirping and bright soliton-higher or- 
der chirping. A mutual cancellation will occur at some 
point (s) when both kinematic chirping and higher order 
chirping are comparable and have opposite signs, and 
will result in chirp reversal at the point (s) of cancella- 
tion. Chirp reversal plays a significant role in fiber op- 
tics, and has attracted considerable attention recently in 
context of pulse retrieval in disp ersion- nonlinearity man- 
aged optical fibres [H, [13, El 
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FIG. 2: Phase profile of dark soliton plotted in Figure 1. (in 
red color). Blue curve shows the contribution from linear 
chirp whereas green curve shows the contribution from higher 
order chirp. The chirp reversal is clearly seen as peaks in the 
red curve. 



against ^ in figure ([2]) where the chirp reversal is clearly 
seen as two maxima in the profile. 

It should be noted that equation ([8]) with 6*4 7^ has 
no kink solutions, which arc of the type: 



A + Btanh(^) 
C + Dt£inh{C) ■ 



(10) 



However, for 04 — 1/4, 64 is zero allowing exis- 
tence of these kind of solutions. In this case, k — 

(91^+292 
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B 
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TD whereT^-M^eS^^' C = 1 and 



D = ± 



2fc-492A-393A2 



2e-i{A-r) — • Figure (O shows intensity pro- 
file of a typical kink solution. These solutions being 
asymmetric around ^ = 0, have interesting phase pro- 
file which show chirp reversal only once along the profile. 
The condition for existence of this reversal is given by: 
, > and -1 < < 1, and the point of 

304+205 — Duc — B — ' ^ 

reversal is £,r = tanh^^( ^^°^^ ). 

Solutions to equation ([5]) via (O are not restricted to 
localized ones alone, periodic solutions also exist. Infact 
it is easy to show that: 
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is a periodic solution of equation 

9fl 2/9/1 _ 993f, , 993!.\ i_ 894('-3«2 93 , 
ZPi ^Z(72 "'" 1694''^ 3 
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Apart from the solutions discussed above, amplitude 
equation ([7]) albeit with different parameters, has been 
carefully studied in context of cubic quintic nonlinear 
Schrodinger equation [2l|. It has been shown that 
this equation possesses a rich solution space, where the 
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solutions are expressible in terms of Weierstrass func- 
tions, and nature of the solution crucially depends upon 
initial conditions. Similar analysis for this system would 
be relevant, and will shed light on the structure of solu- 
tion space. 

A natural question arises whether the model is intc- 
grable in this regime or not. Following the Ablowitz- 
Ramani-Segur algorithm, we investigate singularity 
structure of the ordinary differential equation ([8]) , which 
is obtained from an exact reduction of the original partial 
differential equation ^ [13, HI] . Interestingly, we found 
that ordinary differential equation represented by ([5]) 
possess poles as the only movable singularities, which im- 
plies that this system indeed has the P-property [H, HI] ■ 
Hence, we see that this system passes Painleve test, and 
is Painleve integrable, which guarantees existence of reg- 
ular solutions in general. 

In conclusion, we have found a new class of exact solu- 
tions to NLSE system in the presence of self steepening 
and self frequency shift terms. These include localized 
solutions of dark-bright type, kink solutions and periodic 
solutions. These solutions have nontrivial phase chirping 
which varies as a function of intensity and are different 
from that in Ref. where the solutions had a trivial 
phase. A nontrivial connection of this sytem with linear 
Schrodinger equation in appropriate limits is pointed out. 
A duality is seen between the dark-bright type of solu- 
tion and kinematic-higher order chirping. A novel form 



of self phase modulation has been observed in this case, 
which shows chirp reversal across the pulse profile. It 
is known that prechirping of pulses often leads to better 
quality of pulses, in particular it is quite effective with 
distributed GVD and nonlinearity [H, [H, [H]. In 
this context, the solutions having chirping due to initial 
conditions as well as dynamical conditions will provide a 
better control. It is noted that for some parameter val- 
ues the intensity and phase of these solitons will exactly 
be the same as NLSE solitons, and are found to be chi- 
ral, with the direction of propagation controlled by self 
steepening term. Both dark-bright solitons are found to 
exist in normal- anomalous dispersion regimes. It is seen 
in some cases that intensity of these solitons would be 
like NLSE soliton and only phase structure will be dif- 
ferent. Kink solutions are found to exist in this system 
for special choice of parameters. The system is seen to 
possess P-property and hence is Painleve integrable. 
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